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We report on a search for CP violation in the decays D ± —> K®K ± , 
Df — >• K^K ± and Df — >• using a data set corresponding to an 

integrated luminosity of 469 fb _1 collected with the BABAR detector at 
the PEP-II asymmetric energy e + e~ storage rings. The CP- violating 
decay rate asymmetries A C p are determined to be (+0.13 ± 0.36(stat) ± 
0.25(syst))%, (-0.05±0.23(stat)±0.24(syst))%, and (+0.6 ± 2.0(stat) ± 
0.3(syst))%, respectively for the three modes, before correction for the 
relevant asymmetry. After this correction, all the measurements are 
consistent with zero within one standard deviation. These are currently 
the most precise measurements of these asymmetries. 
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1 Introduction 



CP violation (CPV) in charm decay is a sensitive probe of physics beyond the Stan- 
dard Model (SM). Owing to its suppression within the SM, a significant observation of 
direct CPV in charm decay would indicate the possible presence of new physics effects 
in the decay process. In a previous publication PQ, we reported a precise measurement 
of CP asymmetry in the decay D ± —> K®^, where the measured asymmetry was 
found to be consistent with the value expected for indirect CPV in the K° system. 
The LHCb and CDF collaborations have recently reported evidence for CPV in the 
decays D° —> K + K~ and ,0° —> 7r + 7r _ [2 [3], which, if confirmed, could be due either 
to new physics effects or to significant enhancement of penguin diagrams in charm 
decays [H [5] . Additional information and corroboration of these observations in other 
channels are necessary to resolve these questions. 

In this report we present measurements of direct CPV asymmetries in the decay 
channels: D ± -> K^K*, Df -> K^K*, and Df -> K^. For these channels, 
as in the case of the decay D ± —> K®^, we expect a K°-induced asymmetry of 
« (±0.332±0.006)% [6]. The sign of the X°-induced asymmetry is positive (negative) 
if a K° (K°) is present in the corresponding tree level Feynman diagram. The exact 
magnitude of the asymmetry would depend on the requirements on the reconstructed 
7r + 7r - decays and the decay kinematics [7j. Previous measurements of Aqp in 
these channels have been reported by the CLEO-c [8] and Belle collaborations j9]. 

For this analysis we employ a technique similar to that used in the measurement 
of CPV in D ± — > [I], thus for the description of certain analysis details we 

refer to our previous publication. 

2 Data Selection and Fit Procedure 

The data used for these measurements were recorded at or near the T(AS) resonance 
by the BABAR detector at the PEP-II storage rings. The BABAR detector and the 
coordinate system used throughout are described in detail in Refs. jTQl [TT]. The 
total integrated luminosity used is 469 fb _1 . To avoid any potential biases in the 
measurements, we adopt a "blind analysis" approach, where for each channel we 
finalize the whole analysis procedure prior to extracting Acp from the data. 

Signal candidates were reconstructed by combining a candidate reconstructed 
in the decay mode 7r + 7r~ with a charged pion or kaon candidate. A candi- 

date is reconstructed from two oppositely charged tracks with assigned the charged- 
pion mass. candidates are required to have an invariant mass within ±10 MeV/c 2 
of the nominal mass [6], which is equivalent to slightly more than ±2.5 a in the 
measured mass resolution. The x 2 probability of the tt + tt~ vertex fit must be 
greater than 0.1%. To reduce combinatorial background, we require the measured 
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flight length of the candidate to be greater than 3 times its uncertainty. A re- 
constructed charged-particle track that has pt > 400 MeV/ c is selected as a pion or 
kaon candidate, where pt is the magnitude of the momentum in the plane perpendic- 
ular to the e + e~ collision axis. At BABAR, charged hadron identification is achieved 
through measurements of specific ionization energy loss in the tracking system, and 
the Cherenkov angle obtained from a detector of internally reflected Cherenkov light. 
A CsI(Tl) electromagnetic calorimeter provides photon detection, electron identifica- 
tion, and neutral pion reconstruction [10]. In our measurement, it is required that 
a pion candidate not be identified as a kaon, a proton, or an electron, while a kaon 
candidate is required to be identified as a kaon and not as a pion, a proton, or an elec- 
tron. The criteria used to select pion or kaon candidates are very effective in reducing 
the charge asymmetry from track reconstruction and identification, as inferred from 
the study of control samples, as described below. A kinematic vertex fit to the whole 
decay chain is then performed, and it is required that each candidate form a vertex 
close to the beam interaction region [12J. We retain only candidates having a x 2 
probability for this fit greater than 0.1%, and an invariant mass m(K®h), h = 7r, K ) 
within ±65 MeV/ c 2 of the nominal mass [6j; this is equivalent to more than ±8 a 
in the measured Df, mass resolution. 

We further require the magnitude of the Df candidate momentum in the e + e~ 
center-of-mass (CM) system, p*(Df), to be between 2.6 and 5.0 GeV/c to suppress 
the combinatorial background from BB events. For the decay channel D ± — >> K®K ± ) 
which suffers from lower statistics, we retain candidates with p*(D ± ) between 2.0 and 
5.0 GeV/c. This selects some signal events originating from B meson decays [13], while 
maintaining an acceptable level of combinatorial background. Additional background 
rejection is obtained by requiring that the impact parameter of the candidate with 
respect to the e + e~ interaction region [TO] , projected onto the plane perpendicular 
to the collision axis, be less than 0.3 cm, and that the lifetime r xy (D^) be 
between —15 and 35 ps. The lifetime is measured using L xy (Df\ defined as the 
distance of the decay vertex from the interaction region projected onto the plane 
perpendicular to the collision axis. 

In order to further optimize the sensitivity of the Acp measurements, a multi- 
variate algorithm is constructed from seven discriminating variables for each 
candidate. These are: T xy (D^), L xy (D^), p*(D^), the momentum magnitude and 
transverse component with respect to the beam axis for the and for the pion 
or kaon candidate. For D ± — >> K®K ± and Df — >> K®K ± the multivariate algorithm 
with the best performance is a Boosted Decision Tree (BDT), while for Df — » 
the best algorithm is a Projective Likelihood (LH) method [14]. The final selection 
criteria, based on the output of the relevant multivariate selector, are optimized using 
the S/^/S + B ratio as defined in Ref. pQ. 

For each mode the signal yield is extracted using a binned maximum likelihood 
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Figure 1: Invariant mass distribution for (a) D ± —> K®K ± , (b) Df K®K ± ) and 
(c) Df Kgiv^ candidates from the data (black points). In each figure, the solid blue 
curve shows the fit to the data, the dashed red curve is the sum of all backgrounds, 
while the dotted green curve represents combinatorial background only. 

(ML) fit to the distribution of invariant mass, m^K^h), for the selected candi- 
dates. The total probability density function (PDF) is the sum of signal and back- 
ground components. The signal PDF is modeled as a sum of two Gaussian functions 
for both the D ± — >> K®K ± and Df —> K®K ± modes, and as a single Gaussian function 
for Df —> Kgiv^ mode. The background PDF is taken as the sum of two components: 
background candidates from misidentified and partially reconstructed charm decays, 
and a combinatorial background component from other sources. For D ± K®K ± 
(Df K®^) the charm background is mainly from the tail of the invariant mass 
distribution of Df (D ± K^) candidates. For Df the 

charm background originates from D ± —> where the ^ is misidentified as 

K ± . This wrong mass assignment shifts the reconstructed invariant mass to higher 
values, and the resulting distribution is a broad peak with mean close to the Df 
mass. The charm background distributions are modeled using a PDF sampled from 
the simulation events for these modes. The combinatorial background is described as 
a first (second)-order polynomial for Df (D ± -> K^K* 1 and Df -> K^K^. 

The fits to the m(K° s h) distributions yield (159.4 ± 0.8) x 10 3 D ± -> K^K* decays, 
(288.2±1.1) x 10 3 Df decays, and (14.33±0.31) x 10 3 Df decays. 

The data and the fits are shown in Fig. [TJ All of the fit parameters are extracted 
from the fits to the data. 



3 Aqp Asymmetry Extraction 

For each channel, we determine Aqp by measuring the signal yield asymmetry A 
defined as: 

Nd+ ~ N D - 

A = N l \ N W > (!) 
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We consider A to be the result of two other contributions in addition to Acp, namely 
the forward-backward (FB) asymmetry (A FB ) and a detector-induced component. 
We measure A FB together with Aqp using the selected dataset, while we correct the 
data for the detector-induced component using coefficients from a control sample. 

In this analysis we use a data-driven method, described in detail in Ref. [I], to 
determine the charge asymmetry in track reconstruction as a function of the magni- 
tude of the track momentum and its polar angle. The method exploits the fact that 
the T(4aS) — >■ BB events provide sample of evenly populated positive and negative 
tracks, free of any physics-induced asymmetries, hence allowing the determination of 
detector-related asymmetries in the reconstruction of charged-particle tracks. Start- 
ing from a sample of 50.6 fb _1 of data collected at the T(45) resonance and an 
off-resonance data set of 44.8 fb _1 , we obtain a large sample of charged tracks and 
apply the same charged pion or kaon track selection criteria used in the reconstruction 
of K®K ± and Df decays. Then, by subtracting the off-resonance 

sample from the on-resonance sample, we obtain a sample of more than 120 mil- 
lion pion candidates and 40 million kaon candidates, originating from T(4S f ) decays. 
These are then used to compute the efficiency ratios for positive and negative kaons 
and pions. The ratio values and associated statistical errors are shown in Fig. [2j For 
K®K~ (D~ —> K®tt~) decays, the (D~) yields, in intervals of kaon (pion) 
momentum and cos#, are weighted with the kaon (pion) efficiency ratios to correct 
for the detection efficiency differences between K + and K~ (7r + and 7r~), leaving only 
FB and CP asymmetries. The largest correction is around 1% for pions and 2% for 
kaons. 

Neglecting the second-order terms that contain the product of Aqp and A FB) the 
resulting asymmetry can be expressed simply as the sum of the two. Given that 
A FB is an odd function of cos#^, where 6* D is the polar angle of the candidate 
momentum in the e + e~ CM frame, Aqp and A FB can be written as a function of 
I cos#nl as follows: 



where A{+\ cos#^|) is the measured asymmetry for the candidates in a positive 
cos 6^ bin and A(— \ cos6^|) in its negative counterpart. 

A simultaneous Maximum Likelihood (ML) fit to the and D^, invariant mass 
distributions is carried out to extract the signal yield asymmetries for each one of the 



A FB (\cos0* D \) 



A(+\cos0* D \) - A(-\cos0* D \) 
2 



(2) 



and 



A CP (|cos^|) 



A(+|cos^|) + A(-|cos^|) 
2 



(3) 
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Figure 2: Ratio between the detection efficiency for positive and negative tracks 
(top) and corresponding statistical errors (bottom). Plots on the left are for tt + and 
7T _ , while on the right for K~ and K + . The values are computed using the numbers 
of positive and negative tracks in the selected control sample. 



ten equal intervals of cos#^, starting with interval defined as —1.0 > cos#^ > —0.8. 
The PDF shape that describes the distribution in each sub-sample is the same as 
that used in the fit to the full sample, but the following parameters are allowed 
to float separately in each sub-sample (referred to below as split parameters): the 
yields for signal, charm background and combinatorial candidates; the asymmetries 
for signal and combinatorial candidates; the width and the fraction of the Gaussian 
function with largest contribution to the signal PDF; the first order coefficient for the 
polynomial of the combinatorial background. For D ± — >> K®K ± the yield of charm 
background candidates has been fixed to in intervals 0, 1, and 2 in order to obtain a 
fully converging fit. Interval 9 has the lowest number of candidates compared to the 
other intervals, so, for D ± —> K®K ± and Df —> K®K ± ) we use only one Gaussian 
function in this interval for the signal PDF by fixing the fraction of the first Gaussian 
function to 1. For the CPV asymmetry of charm background candidates we use the 
same floating parameters as for the signal candidates, because the largest source of 
CPV asymmetry for both samples is the CPV contribution from K^—K mixing. For 
the mode Df —> i^vT^, where the primary charm background channel, D ± —> i^vT^, 
has the same magnitude but opposite sign asymmetry from K° — K° mixing, we use 
a separate parameter for the asymmetry of the charm background candidates. If 
the fit values for a split parameter are statistically compatible between two or more 
sub-samples, the parameter is forced to have the same floating value among those 
sub-samples only, in order to achieve a more stable fit. For Df ->■ K°tt ± the width 
of the first Gaussian function for the signal PDF is set to the same floating value 
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Figure 3: A CP asymmetry for (a) D ± -> K^K*, (b) Df -> K®^, and (c) D 
K^n ± as a function of |cos6^| in the data sample. The solid line represents the 
central value of Acp and the gray region is the ±la interval, both obtained from 
a x 2 minimization assuming no dependence on |cos#^|. The corresponding A FB 
asymmetries are shown in (d), (e), and (f). 



for intervals 0, 1, 2, and 3. The first order coefficient of the polynomial describing 
the combinatorial background is set to the same floating value for intervals 3 to 
7 (D ± -> K^K*), for intervals 4 to 7 (Df -> K^K*), and for intervals 1 to 6 
(Df -> K^). The final fits for D ± -+ K^K*, Df -> K^K*, and Df -> 
involve a total of 70, 80 and 64 free parameters, respectively. 

The Acp and A FB values are shown in Fig. |3j The weighted average of the five 
A CP values is found to be: (0.16 ± 0.36)% for D ± -> K^K*, (0.00 ± 0.23)% for 
Df — >> K®K ± ) and (0.6 ± 2.0)% for Df — >> K {) s ti ± ) where the errors are statistical 
only. 



4 Analysis Validation and Systematics 

We perform two tests to validate the analysis procedure for each channel. The first 
involves generating 5000 toy MC experiments using the PDF obtained from the fit to 
data, and then extracting Acp from each experiment. For the modes D ± —> K®K ± 
and Df K^K*, we find small biases of -0.036 ± 0.014 and +0.041 ± 0.014 (in 
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Table 1: Summary of the systematic error contributions to the Aqp measurement in 
each mode. The values are given as absolute percentage (BABAR Preliminary). 



Systematic uncertainty [%] 


D ± -> IQK* 


Df K^K* 




Efficiency of PID selectors 


0.05% 


0.05% 


0.05% 


Statistics of the control sample 


0.23% 


0.23% 


0.06% 


Mis-identified tracks in the control sample 


0.01% 


0.01% 


0.01% 


Binning in cos © 


0.04% 


0.02% 


0.27% 


K°—K° regeneration 


0.05% 


0.05% 


0.06% 


K®— K® interference 


0.015% 


0.014% 


0.008% 


Total 


0.25% 


0.24% 


0.29% 



units of statistical error), in the fitted values of the Acp parameter. Therefore we 
apply a systematic correction of +0.013% to the value of Acp for D ± —> K®K ± and 
of -0.01% to that for Df —> K®K ± , to account for this effect. The fit returns an 
accurate estimate of the statistical uncertainty for all the modes. The second test 
involves fitting a large number of MC events from the full BABAR detector simulation. 
We measure Acp from this MC sample to be within ±1<t of the generated value of 
zero. 

A list of systematic error contributions and the quadrature sum for each mode 
is reported in Table [T] Primary sources of systematic uncertainty are the statistical 
uncertainties in the detection efficiency ratios used to weight the yields and the 
contributions from misidentified particles in the data control sample used to determine 
the charge asymmetry in track reconstruction efficiency. 

The technique used here to remove the charge asymmetry from detector-induced 
effects produces a small systematic uncertainty in the measurement of Acp due to the 
statistical error on the relative efficiency estimation. We perform simultaneous fits to 
500 samples of selected Df — >> K®K ± candidates where the applied corrections are 
smeared according to their errors, to produce a distribution of the deviations of Acp 
values from the nominal one. The standard deviation of this distribution is 0.23%, 
and it is assigned as systematic uncertainty to the Df —± K®K ± and D ± — >> K$K ± 
modes, because this contribution depends only on the type of charged particle in the 
mode itself. For the Df — )► mode we assign a systematic error contribution 

obtained as in Ref. pp. This is the dominant source of systematic uncertainty for the 
D^ —> K®K ± modes, as shown in Table jl, 

The small fraction of mis-identified particles in the generic track sample can in- 
troduce small biases in the estimation of the efficiencies, and subsequently in the Acp 
measurements. Because of the good agreement between data and MC samples, we can 
use the MC-simulated candidates to measure the shift in the Aqp value from the fit 
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when the corrections are, or are not, applied. Again, this contribution depends only 
on the type of charged-particle track, hence for the Df —> mode we assume 

the same shift as that obtained in Ref. PQ, namely +0.05%. Fitting the Df —> K®K ± 
MC sample when the corrections are, or are not, applied, we obtain a shift of +0.05%, 
and we assume this same uncertainty for the Df K®K ± and D ± —> K®K ± modes. 
As a result, for each mode we shift the measured Aqp value by the net systematic 
shift to correct for this bias, and then, conservatively, include the same value as a 
contribution to the systematic uncertainty. 

Using MC simulation, we evaluate an additional systematic uncertainty of ±0.01% 
due to a possible charge asymmetry present in the control sample before applying the 
selection criteria. Another systematic uncertainty from the simultaneous ML fit is 
due to the choice of interval-size in cos#^. This can be estimated using the largest 
Acp deviation when the fit is performed using 8 or 12 intervals of | cos#^|, instead of 
10. This is the dominant source of systematic uncertainty for the Df —> mode, 
as shown in Table [U 

We also consider a possible systematic uncertainty due to the regeneration of 
neutral kaons in the material of the detector, since K° and K° mesons produced 
in the decay process can interact with the material in the tracking volume before 
they decay. Following a method similar to that described in Ref. [15], we compute 
the probability for a K° or K° to interact inside the BABAR tracking system, and 
estimate an associated systematic uncertainty of 0.05 — 0.06%. 

Another systematic contribution can be generated by the interference between the 
amplitudes of intermediate and K®, because we are reconstructing our using 

—> 7T + 7r _ , which is not a pure amplitude [7j. The bias in the asymmetry contri- 
bution induced by K°—K° mixing depends directly on the reconstruction efficiency 
as a function of the proper time and, approximately, the more the efficiency is not 
constant, the larger becomes the bias. We produce the K® reconstruction efficiency 
distribution as a function of the proper time using MC truth-matched events after 
the full selection. Then, following the method in Ref. [7], we estimate the asymmetry 
correction factor AAqp defined as: 

AAcp = A C SF - A%>, (4) 

where A^ p is the value obtained from the fit and A^p is the corrected value. The 
correction factors are reported in Table [2] and, to be conservative, we include their 
absolute values as a contribution to the systematic uncertainty. In a similar mode, we 
also estimate the correction factor for D ± —> K^tv^ mode using the reconstruction 
efficiency distribution after the selection detailed in Ref. [lj and we find a value of 
+0.002%. All these corrections are rather small, even compared to those estimated 
in a similar analysis [T6] . The smaller values of the corrections in this analysis are 
due to the improved efficiency for mesons with short decay times that we obtain 
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Table 2: Summary table for Acp measurements. Uncertainties, where reported, are 
first statistical, and second systematic (BABAR Preliminary). 





D ± -> IQK* 


Df -+ K^ 


Df -+ K"** 


Acp value from the fit 


(0.16 ±0.36)% 


(0.00 + 0.23)% 


(0.6 + 2.0)% 


Bias Corrections 




Toy MC experiments 

PID selectors 

K®— K® interference 


+0.013% 
-0.05% 
+0.015% 


-0.01% 
-0.05% 
+0.014% 


-0.05% 
-0.008% 


Acp corrected value 


(0.13 + 0.36 + 0.25)% 


(-0.05 + 0.23 + 0.24)% 


(0.6 + 2.0 + 0.3)% 


Acp contribution 
from K° — K° mixing 


(-0.332 + 0.006)% 


(-0.332 + 0.006)% 


(0.332 + 0.006)% 


Acp value (charm only) 


(0.46 + 0.36 + 0.25)% 


(0.28 + 0.23 + 0.24)% 


(0.3 + 2.0 + 0.3)% 



by making a requirement on the decay length divided by its uncertainty rather than 
on the decay length alone. 

5 Conclusions and Acknowledgements 

In conclusion, we measure the direct CP asymmetry, Acp, in D ± K®K ± , Df 
K^K*, and Df decays using approximately 159,000, 288,000, and 14,000 

signal candidates, respectively. The measured Acp value for each mode is reported 
in Table [2j where the first error is statistical and the second is systematic. In the last 
row of the table, we also report the Aqp values after subtracting the expected Acp 
contribution in each mode due to K° — K° mixing. These results are consistent with 
zero and with the SM prediction within one standard deviation. 
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